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Proper times associated with Milky Way's rotation curve have been computed outside the frame- 
work of general relativity. It is proposed that the flat Minkowski metric can deviate in different ways 
for different two body systems. This effect is more pronounced on galactic scale due to large vari- 
ations in gravitational potential caused by non-uniform distribution of galactic matter. When the 
proper times of stars which are functions of the deviation factors, and the virial mass of the galaxy 
are introduced in modified Newtonian dynamics, they yield precise values of circular velocities for 
the stars of the Milky Way. So the formalism becomes comparable to MOND formalism except that 
the virial mass is used and the MOND constant acceleration parameter gets replaced by a variable 
expression related to the proper time. Hence it becomes possible to extract an equivalent variable 
acceleration parameter which is of the same order of magnitude as the constant MOND parameter. 
This leads to the conclusion that the MOND theory provides an ad hoc empirical expression for the 
effect of proper time on Newtonian dynamics. The expression for proper time derived here when 
introduced in Kepler's third law gives us a modified law that satisfies both the solar system dynam- 
ics as well as the galactic rotation curves. The discussion is focused on Milky Way but the theory 
can be extended to other galaxies. Gravitational redshift of light, bending of light and perihelic 
precession of planets are within the permissible limits and this part of the theory called the periodic 
relativity (PR) also satisfy Einstein's field equations. The theory yields a very accurate and precise 
value for the anomalous acceleration of the Pioneer spacecraft. 
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1. INTRODUCTION 

In Solar System, the planets at the outer edge of 
the system have a lower orbital velocity than those closer 
to the Sun. When the observations of the rotation rates 
of the spiral galaxies began in late 70s, it became clear 
that the galaxies were not rotating in the same manner 
as the Solar System. Stars near the outer edge of the 
galaxy were orbiting at the same speed as stars near the 
middle. Astronomers described this rotation pattern as 
the "flattening of galaxies' rotation curves". Two group 
of theories emerged to explain the flattening of rotation 
curves. One group explained the phenomenon as due 
to the presence of unseen matter or cold dark matter 
(CDM) [E Q, __J and the other group proposed modifi- 
cation to Newtonian dynamics at very large distances 

Some of the related literature on the subject in- 
cludes following. CDM articles descri bing universal ro- 
tation curves 0, F(R)-gravity mode l 1431. b rane world 
model H3, dark matter distribution [TElII QJ, Q3, EE 
E_l , scalar-tensor theories [l3[ , gravitational suppression 
model [14] , scaling relations [la ]. Milky Way rotation 
curves [3 [H, f2p|. |28|. No dark matter articles describ- 
ingMOND theories]!, El El, JM El and other theories 
@rtl, H, H3, |32L Data articles (ml l3 EI M, EE El and 
test articles [H, El EE EE ElUl • 

The theory presented in following sections is based 
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on earlier work [371 ] where I found it necessary to rede- 
fine the orbital energy of the body and eliminate weak- 
field approximation from the solution to Einstein's field 
equations. In this theory the deviation introduced to 
flat Minkowski metric (and therefore the proper time of 
the body) is a free parameter, which can be locked only 
by experimental measurement of proper time, direct or 
indirect. Therefore every two body system can deviate 
differently unlike the general relativity and the proper 
time of the body becomes a critical parameter like mass, 
velocity, radial distance etc. When this deviation is small 
enough as in case of gravitational redshift of light, bend- 
ing of light and the perihelic precession of planets, the 
theory can be easily written in a metric form to the first 
order accuracy. In case of Pioneer anomaly [38| the de- 
viation is much larger but within the range of metric 
representation. In case of galactic rotation curves these 
deviations are so large that it is not possible to write 
this theory in metric form without introducing large er- 
rors. This is because the higher order terms cannot be 
ignored. Therefore the mathematical treatment of the 
rotation curves part of this theory is more or less like the 
special relativity. 

The theory presented in this paper simultaneously 
proposes cold dark matter with modified Newtonian dy- 
namics which transits smoothly from short range to long 
range without abrupt discontinuity and all the parame- 
ters of the theory have a physical meaning. In addition, 
the smallest possible ups and down in what is generally 
described as the flat rotation curve can be accounted by 
the theory. The theory can be written in CDM format 
or MOND format. 
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2. PHYSICS BEHIND THE PRESENT THEORY 

In earlier work named the periodic relativity (PR) 
[37l ] it was proposed that the time is a periodic phe- 
nomenon and can only be defined as period between two 
events. The most fundamental aspect of time is the 
period of a wave. We defined a general form of Max 
Planck's quantum hypothesis applicable to massive par- 
ticles and derived a relation that showed that when a 
force acts on a fundamental massive particle, it acceler- 
ates the particle and simultaneously alters its associated 
de Broglie frequency and therefore force consists of two 
components, the Lorentz force and what we called the 
de Broglie force responsible for effecting the change in 
the associated frequency of massive particle [13 ■ When 
de Broglie force is ignored, the invariant relationship be- 
tween force and energy is destroyed. For photon tra- 
jectory in a gravitational field, we made use of this de 
Broglie force component to derive the gravitational red- 
shift of light without making use of the weak field ap- 
proximation and we made use of the Lorentz force com- 
ponent to derive the deflection of light, again without 
making use of the weak field approximation. In case of 
a gross body such as a planet, we can convinicntly deal 
with the Lorentz force but it is not possible to directly 
account for the de Broglie force. To compensate for this 
deficiency, we make use of the concept of proper time in 
PR. Hence PR [13, proposes a definite connection be- 
tween the proper time of a gross body and the associated 
de Broglie wavelengths and frequencies of its constituent 
fundamental particles which are in bound states. 
We have the fundamental relations E = mc 2 = hv. 



m Vq dr ^1 - /3 2 ' 

\ dv 1 dm d / dt \ dj v dv ^ 
vq dt mo dt dt \dr J dt c 2 dt 

Therefore above relations show that the change in as- 
sociated de Broglie frequency is same as the change in 
relativistic mass of the particle. Altered energy level of a 
body would imply altered relativistic mass of the body. 
Since in PR, the gravitational attraction exists between 
the relativistic masses, the attractive force between the 
two bodies would also be altered. The third term in 
Eq. (|2.2|) would imply rate of change in the ratio of the 
coordinate time and the proper time. The fourth term 
would imply a rate of change in the flat Minkowski met- 
ric and the last term relates all above changes to the rate 
of change of the velocity. 

Single most major difference between PR and the 
general relativity is that the later allows only one time 
specific deviation to the flat Minkowski metric which we 
call the weak field approximation, where as PR allows 
different deviations to the flat Minkowski metric to suit 
different two body problems. Therefore there is no such 
thing as a spacetime fabric or rubber in PR [39| . 



Now how do we apply these concepts to the prob- 
lem of 'flat rotation curves'. First we need to imagine 
how many trillions and trillions of fundamental parti- 
cles have gone into making of a star. The gravitational 
potential within the virial radius of a galaxy does not 
change as uniformly and smoothly as it would in case of a 
spherically symmetric body such as a star. Every galac- 
tic circular radius is associated with a specific gravita- 
tional potential which is dependent on the galactic mass 
distribtion. The state of the associated de Broglie fre- 
quencies (or wavelengths) of the individual constituent 
massive particles (which are in bound state) of the star 
arc directly dependent on the local gravitational poten- 
tial. This state of the constituent particles of the star 
defines the deviation to the flat Minkowski metric with 
respect to the two body system consisting of the center 
of glactic mass and the star. This in turn fixes a proper 
time for the star at the given circular radius with respect 
to the coordinate time of the galactic center. For a par- 
ticular star in circular orbit there is no frequency shift 
of constituent particles and therefore the orbital energy 
is constant but as we move radially to other orbits and 
other stars we find that they have fixed their proper time 
in the same manner depending on the local gravitational 
potential. The energy levels of the galactic stars depend 
on the associated frequencies of their constituent parti- 
cles which define the relativistic mass of the star which 
through Lorentz transformation define the velocity of the 
star which is dependent on the proper time of the star. 
Earlier we saw that the gravitational force consists of 
two components. The Lorentz force and the de Broglie 
force. In short range galacto dynamics, the Lorentz force 
dominates and the de Broglie force is negligible. At long 
range the Lorentz force is negligible and the de Broglie 
force dominates. This statement when analyzed in terms 
of proper time will show us in following section that New- 
ton's second law gets altered by a factor equal to square 
of the ratio of proper time to coordinate time. 



3. EQUATIONS OF STELLAR MOTION 

A. CDM format 

Following is the extention of the theory developed 
earlier [U where we introduced deviation to the flat 
Minkowski metric due to the gravitational field in the 
form, 

(^) 2 = 7 2n = (i-/? 2 r\ (3.1) 

For small values of (3 and n, Eq. (|3.1[) can be written in 
approximate form to first order accuracy as 

1 _ 2^J ' (3 - 2) 
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where t is the coordinate time, r the proper time of the 
orbiting body and n is a real number. The corresponding 
line element in polar coordinates is, 



ds 2 



- 2 dt 2 



ndr' 



nr 2 d6 2 



n(r 2 sin 2 0)d<f> 2 . 



(3.3) 



We showed that the line clement Eq. (|3.3|) satisfies Ein- 
stein's field equations for any constant value of n or for 
n = 0. However, for very large values of n, approximate 
relations given by Eqs. (|3.2[) and (|3.3[) are invalid. In pre- 
vious analysis, we have used n = for photon trajectory, 
n = 1 for flat Minkowski metric, n = 4 for solar system 
planets. These values are small enough to allow the met- 
ric form Eq. (|3.3[) . In case of Pioneer anomaly, n is of 
order 10 5 which is reasonably small when [3 2 is of order 
10 -9 and metric representation is permissible. In case of 
Milky Way rotation curves as we will see, n is of order 
10 6 and 1 is of order 10" 7 and if Eqs. (HHJ) and ([Ojl 
are used, they would introduce large errors. Therefore 
we will restrict our use to Eq. (|3.1[) . Under the circum- 
stances we also conclude that the theory developed below 
cannot satisfy Einstein's field equations even if n is a con- 
stant of circular orbit. 

In our earlier analysis [37] ] of perihelic precession of 
planets, we developed the modified Newtonian dynamics 
expression 



cPr 
dt 2 



(3.4) 



We can use Eq. (|3.1[) to introduce proper time in 
Eq. (I3~D . 



d 2 r _ pL f v 2 
dr 2 r 



r = 



(3.5) 



Here cj>(t) is the gravitational potential w.r.t. the coordi- 
nate system of the galaxy and 4>(t) is the gravitational 
potential w.r.t. the coordinate system of the star. 

Since the transverse component of the gravitational 
acceleration is absent in the polar coordinate system of 
the central potential, we can write 



dr 2 



d 2 r 
dr 2 



dr 2 r 3 J r 2 \ 

For circular orbits we have 

and 



(l-n) 



(3.6) 



(3.7) 



dr 2 



h = 



Substitution of Eq. (|3~8|) in Eq. J3J|) gives 



v 
r 



(l-n) 



(3. 



(3.9) 



Eq. (|3.9p can be written in the logarithmic form as 



= In 1 - 



— n In I 1 — 



In 



From Eq. (|3.10[) we get the expression for 



(3.10) 



1 - 



ln( S-L 



lnfl 



(3.11) 



Substitution of Eq. (|3.11|) in Eq. (|3.ip gives the expres- 
sion for proper time 



dr = dt I 1 - — 



1 12 



(3.12) 



When virial mass and circular velocities are introduced 
in Eq. (|3.11[) it gives deviation factor n and proper time 
interval dr which corresponds to the CDM theory. 



1. CDM MOND connection 

In order to relate the CDM theory of sec.3A with 
the MOND formalism we can write Eq. (|3 . 9[) in the form 



v 
r 



GM D (M v 
r 2 \Md 



2\ (l-n) 



(3.13) 



where M v is the virial mass (halo mass) and is the 
stellar mass (disk mass). 



f_ ( M v 
r \M D 



2\ (I"") 



1 - 



GMt 



(3.14) 



The L.H.S. of Eq. (j3~T4|) is equivalent to the MOND ex- 
pression and therefore we can write 



r \M D 



2\ (1-™)' 



e(-)a. (3.15) 



I — ) a = 

' \ao) a 0v /l + (a/a ) 2 ' 



v 
r 



My 

M D 



1 - 



(l-n) 



a 
a- 



(3.16) 



(3.17) 



Therefore the expr ession a e is equivalent to the MOND 
function a m = ao yl + (a/ao) 2 . This means that if we re- 
place the MOND function by a e , both CDM and MOND 
theories would yield precisely the same rotation curves 
with same proper time intervals. In addition, the arbi- 
trary MOND function will now have a physical meaning. 
Some of the difficulties encountered by MOND in fitting 
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the Milky Way curve near solar radius [23j, |26[ gets elim- 
inated. Same can be said about NGC3198 and clusters 
of galaxies [Z(| • 

With respect to rotation curves of other galaxies, 
if one knows the stellar mass Mb [3, the virial mass 
M v = Md + Mh can be computed using the relation 
given by Shankar et al. [36|. Here Mh is the halo mass. 



M D = 2.3 x 1O-M [M * /(3 X 10llM ° )] 



3.1 



1 + [M h /(3 x 10"M Q )] 2 - 2 
From Eqs. ([3~T7]) and (|3~T2"]) we can write 



(3.18) 



r \ c 2 



{dr/dtf 



V 2S 


\ M v 




l M D 


v 2 \ 


M v 




M D 


M v 





{dr/dt) 2 M L 



TIT- ( 3 - 2 °) 



(3.21) 



Eq. (|3.21j) provides partial justification to the MOND 
claim that Newtons second law gets modified in the deep 
MOND range. But there is no justification for the pres- 
ence of the mass ratio M v /Mry- In the central bulge 
of the galaxy, factor (dt/dr) 2 (M v /MD) seems to remain 
close to unity which would mean that there is no viola- 
tion of Newton's second law. In outer flat disk area, this 
factor deviates from unity implying violation of Newton's 
second law. Looking at it from CDM point of view, there 
is no violation of Newton's second law when the velocity 
and acceleration are computed in proper time of the star 
and there is no need of introducing the ratio M v /Md- 
This would also mean that the line of sight velocities 
of stars measured by the astrophysicists using Doppler 
shifts of their spectral features arc actually the velocities 
measured in proper time of the stars. 



2. Proper time intervals of stars in Milky Way 

In this section we apply the equations of stellar mo- 
tion developed in previous section to Blue Horizontal- 
Branch (BHB) halo stars of the Milky Way [Hj]. The 
circular velocity estimates are based on Naab's simula- 
tion [4lJ for which Xue et al [l^| obtained virial mass 
M v = 1.12 ± 0.17 x 1O 12 M corresponding to the adia- 
batically contracted NFW profile. To this data, one ad- 
ditional data point for solar radius of 8kpc [26| is added 
and the results are shown in Table |TJ The first thing 
we notice is that for all values of r, n is negative imply- 
ing that time runs faster and faster as one moves away 
from the glalactic center till the peak velocity is attained 
on the rotation curve. From here as r increases there is 
slight drop in the velocity which corresponds to the slow- 
ing down of time. Secondly, functions a e and Qui 3XC of 
the same order of magnitude and a e provides the perfect 
fit which a m cannot. 



TABLE I: Milky Way data shows time running faster at solar 
radius by 3.5 times and slows down at larger radii. a e and a m 
in m/s 2 . CDM format. 



r(kpc) 


v(km/s) 


n 


dr/dt 


a e x 10 10 


a m x 10 10 


7.5 


216 


-5051271 


3.71 


3.22 


2.34 


8.0 


220 


-4681268 


3.527 


3.466 


2.29 


12.5 


227 


—3509342 


2.73 


3.93 


1.796 


17.5 


179 


-6032730 


2.93 


1.519 


1.34 


22.5 


168 


-6452233 


2.75 


1.18 


1.267 


27.5 


183 


-4440256 


2.287 


1.66 


1.263 


32.5 


143 


-8705548 


2.69 


0.618 


1.217 


37.5 


170 


-4639096 


2.108 


1.236 


1.225 


42.5 


183 


-3271979 


1.84 


1.66 


1.227 


47.5 


165 


-4341245 


1.93 


1.096 


1.214 


55 


183 


-2580035 


1.61 


1.66 


1.216 






3. Kepler's third 


law 





There is a considerable difference between the 
proper time given by Eq. (|3.12p and that given by the 
general relativity formula for Kcplcrian circular orbits. 



dr(gr) = dt [ 1 



3[i 
2c 2 R 



(3.22) 



If we put n — 3 and v 2 = ji/R in Eq. (|3 . 1 2[) we get 
Eq. (|3.22[) implying that GR would allow only one fixed 
value for n for all kind of two body problems. This is 
what we call the weak field approximation. This is how 
GR locks the proper time in the theory which is based 
on the near earth observations. As long as the proper 
time data for distant objects are not available, no one 
can question the validity of GR. The present theory also 
has the same adavantage. With respect to all Milky Way 
data points given in Table [H ratio dr/dt given by GR is 
approximately of order (1 — 10 -6 ) which is almost unity 
and instead of time running faster with increasing veloc- 
ities, it actually runs slow. 

Coming back to the argument that the line of sight 
velocities of stars measured by the astrophysicists using 
Doppler shifts of their spectral features are actually the 
velocities measured in proper time of the stars. This fact 
also clearly shows up in Eqs. (|3.7[) . (|3.8|) and (|3.9p where 
v = ds/dr. Since Newtonian acceleration is a = (v 2 /r) 
and with such a large ratios for dr/dt as given in Tabled 
at first glance it would appear that Newton's second law 
is violated. If one were to use GR value of proper time 
which is based on Kepler's third law, one would hardly 
notice any effect on Newton's second law. However, it is 
well established that the spiral galaxy rotation curves do 
not fall off as predicted by Kepler's third law. Now Ke- 
pler's third law has to do with the orbital period of the 
body and therefore the time. Hence the study of rotation 
curves basically tells us that our understanding of time 
is not perfect. Towards this goal if we were to rewrite 
Kepler's third law in PR we can make use of Eqs. 



5 



and (|3 . 1 2[) to obtain 
47r 2 7- 2 47r 2 r 2 



P 2 (r) P 2 (t) 



1-^ 



r \ c z 



F 2 fr) 



47r 2 r 3 



P 2 (t) 



47r 2 r 3 



(n-l) 



2^,3 



47r 2 r 



(3.23) 



(3.24) 



(3.25) 



where P(t) and -P(i) are the orbital periods in proper 
time and coordinate time respectively. Therefore the 
data corresponding to the solar radius in Tableland with 
the virial mass of the Milky Way, Eq. ([3~2"I|) yields 223.4 
million years for the orbital period of the sun around the 
galactic center. For n = 4 for the solar system, Eq. (|3.24|) 
yields Kepler's third law as a close approximation. For 
all small values of n, P(t) « P(t). 



TABLE II: Milky Way data with MOND equations show time 
running slower than in CDM theory. MOND format. 



r(fcpc) 


v(km/s) 


n 


dr/dt 


7.5 


216 


610460 


0.853 


8.0 


220 


762648 


0.814 


12.5 


227 


1 365558 


0.676 


17.5 


179 


354588 


0.938 


22.5 


168 


-230011 


1.03 


27.5 


183 


732196 


0.87 


32.5 


143 


-2973816 


1.4 


37.5 


170 


25579 


0.996 


42.5 


183 


810578 


0.86 


47.5 


165 


-336138 


1.05 


55 


183 


834207 


0.85 



of function £(a/ao). The time ratio dr/dt remains closer 
to unity in contrast with the CDM format. Two theo- 
ries present sharp contrast at shorter range. In MOND 
format the abrupt transition from very short range to 
long range cannot be eliminated. In CDM format this 
transition is smooth and continuous. 



MOND format 



CONCLUSION 



If a resourceful theory comes along and shows that 
the dark matter cannot exist and MOND is the best op- 
tion, then if we were to look for a proper time signature 
of MOND, the equations can have following forms. If 
one were to choose Eq. (|3.4p for introducing the MOND 
function Eq. (|3.16[) . the MOND theory would remain as 
it is today and there cannot be any discussion of proper 
time. If one introduces Eq. (|3. 16[) in Eq. (|3.5p instead, 
the equations would be as follows. 



«o\/l + ( a / a o) 



/ 2 \ ( 1— r 



(3.26) 



hi 



v 



r>2 ^,2 



-1/2- 



In 



-1/2' 



ln(l 



(3.27) 



(3.28) 



The proper time expression remains same as Eq. (|3.12[) . 
Values obtained from Eqs. (pHSj) and (gig) for BHB 
halo stars of the Milky Way are shown in Table [Til The 
Milky Way stellar mass used 051] is 5.0924 x 10 10 M Q . 
Value of n remains positive near solar radius and fre- 
quently changes sign mainly due to the empirical nature 



In PR we have proposed that the flat Minkowski 
metric can deviate in different ways for different two body 
systems. The theory is at variance with the formalism of 
proper time used by general relativity. In PR the proper 
time of a body is associated with the gravitational fre- 
quency shift of the constituent fundamental particles of 
the body. The frequency shift changes the energy level 
of the body and alters its relativistic mass which gets 
reflected in its motion. This effect is more pronounced 
on galactic scale due to large variations in gravitational 
potential caused by non-uniform distribution of galactic 
matter. When the proper times of stars which are func- 
tions of the deviation factors, and the virial mass of the 
galaxy are introduced in modified Newtonian dynamics, 
they yield precise values of circular velocities for the stars 
of the Milky Way. So the formalism becomes compara- 
ble to MOND formalism except that the virial mass is 
used and the MOND parameter a m gets replaced by a 
variable expression a e related to the proper time. These 
two parameters are of the same order of magnitude. This 
leads to the conclusion that the MOND theory provides 
an ad hoc empirical expression for the effect of proper 
time on Newtonian dynamics. The expression for proper 
time derived here when introduced in Kepler's third law 
gives us a modified law that satisfies both the solar sys- 
tem dynamics as well as the galactic rotation curves. 

The theory weighs heavily towards the CDM the- 
ory but theories such as MOND cannot be conclusively 
ruled out without direct or indirect experimental verifi- 
cation of the predicted proper times shown in Tables UJ 
fclXTI Gravitational redshift of light, bending of light and 
perihelic precession of planets are within the permissible 
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limits [13] and this part of the theory also satisfy Ein- 
stein's field equations. The theory yields a very accurate 
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